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Abstract: 

We present results for the differential distributions of jets from non-leptonic decays of 
polarized top quarks within the Standard Model, including QCD radiative corrections. 
Our work extends existing results which are only available for semileptonic top quark 
decays at the parton level. For f (f) — > &-jet +2 light jets we compute in particular the 
QCD-corrected top-spin analysing power of the £>-quark jet and the least energetic light 
jet. The dependence of the results on the choice of the jet recombination scheme is found 
to be small. In addition we compute the spin analysing power of the thrust axis. Our 
results constitute a so far missing ingredient to analyse top quark production and subse- 
quent non-leptonic decay at next-to-leading order in a s , keeping the full information on 
the top quark polarization. 
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1 Introduction 



The detailed analysis of the dynamics of top quark production and decay is a major objec- 
tive of experiments at the Tevatron, the LHC, and a possible future linear e + e collider. 
A special feature of the top quark that makes such studies very attractive is its large decay 
width: In contrast to the light quarks the large top decay width « 1.5 GeV serves as a 
cut-off for non-perturbative effects in top quark decays. As a consequence precise theo- 
retical predictions of cross sections and differential distributions involving top quarks are 
possible within the Standard Model and its extensions. A confrontation of such predic- 
tions with forthcoming high-precision data will lead to accurate determinations of Stan- 
dard Model parameters and maybe hints to new phenomena. In particular, observables 
related to the spin of the top quark can be studied and utilized to search for new interac- 
tions of the top quark [[[]]. In e + e collisions, top quarks are produced highly polarized, 
especially if one tunes the polarization of the incoming beams, as possible e.g. at the 
TESLA collider [@]. Furthermore, even for purely QCD-induced production of top quark 
pairs, the spins of t and t are in general highly correlated [J5p. 

The polarization and spin correlations of top quarks must be traced in the differential 
distributions of the decay products. This is possible since the information on the top 
polarization is transferred to the angular distribution of the decay products through its 
weak, parity violating decays. To be more precise: Consider a polarized ensemble of top 
quarks at rest with polarization vector P, < |P| < 1. The differential decay distribution 
with respect to the angle $ between P and the direction p of a given decay product reads: 

1 dT =-(l + lP|K p COSfl). (1) 

Tdcostf 2 y 1 1 p ' v 

In Eq. ([[]), r is the partial width for the corresponding decay of unpolarized top quarks, 
and K p is the so-called spin analysing power of the final state particle or jet under consid- 
eration. For example, in the semileptonic decay t — > l + V[b, the charged lepton (^-quark) 
has spin analysing power K p = +1 (~ —0.41) at the tree level within the Standard Model. 
In hadronic top decays t — > bdu (where d{u) stands generically for J, s (u, c)), the role of 
the charged lepton is played by the d quark. Only for t — ► bsc the maximal spin analysing 
power of the s-quark could in principle be used by tagging b and c-quark jets. However, 
the efficiency of charm-tagging is quite low, and one should try to use a spin analyser that 
is both efficiently detected and has a large analysing power. A good choice is the least 
energetic light (i.e. non-^-quark) jet [Hp, which at tree level has K p = +0.51. This follows 
from the fact that with a probability of 61% (at tree level) this jet contains the J quark. 

The topic of this letter are the QCD corrections to the above tree level results for K p . 
In fact we will be a little more general and discuss corrections to the fully differential 
decay distribution of polarized top quarks to be defined in section ||. These corrections 
are one ingredient for a full analysis of top quark (pair) production and decay at next-to- 
leading order in a s , both at lepton and hadron colliders. They form part of the factorizable 
corrections within the pole approximation [EL 0] for the top quark propagator(s). (For the 
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non-factorizable contributions, see ref. [{7|].) QCD corrections to the production of top 
quark pairs, including the full information about their spins, can be found in ref. pSp for 
collisions and in ref. [Q] for hadron-hadron collisons. In the case of e + e~ collisions, 
also fully analytic results for the top quark polarization [ JTOt ] and a specific spin correlation 
JTT1 ] to order a s are available. In ref. [|l"2"|], helicity amplitudes for e + e — > /fX production 



and semileptonic decays are computed to order a s and are used to construct an event 
generator. The reaction e + e~ — > tt — > W + £V7 _ Z? is treated in ref. [|T3"|], including both 
factorizable and non-factorizable QCD corrections. The theoretical status of polarized top 
quark decay is as follows: A complete calculation of the angular decay distribution for 



f(t) — > to order a 4 can be found in ref. Q14]. The QCD corrections for semileptonic 



polarized top quark decays have been computed in ref. [|T5|]. We will compare our results, 
from which the semileptonic case can be easily derived, to those of ref. [ p3| ] in section |]. 

The outline of this paper is as follows: In the next section we shortly review the tree level 
results for the decay of polarized top quarks. In section [5] we discuss the calculation of the 
QCD corrections. Section [| contains our numerical results, which we discuss in section 

B 



2 Kinematics and tree level results 

Consider an initial state consisting of top quarks at rest with polarization P. For the non- 
leptonic decay 

KPt) HPb)+u(Pu)+d(Pd), ( 2 ) 

the phase space Rj, of the final state may be parametrized by two scaled energies and two 
angles: 

2 

lit 

dRs = 32 ^ 4 dxjdx b d%d cos 0, (3) 

where x^, = 2Eh/m t , xj = 2Ej/m t , cos 9 = P ■ pj, and % is the (signed) angle between 
the plane spanned by P and pj and the plane spanned by pj and p^. We will neglect the 
masses of the light quarks u and d. The differential decay rate is given by 

dT° = —\!M {p b ,p u , Pi )\ 2 dR^ (4) 

where | M (pt, , p u , pj) | 2 stands for the squared matrix element averaged over the colour of 
the initial state and summed over colour and spins of the final state. The fully differential 
distribution for reaction (Q) reads at tree level: 

dT ° ~ ~ ^-^-^ (l + |P|cos0), (5) 



dxjdx b d%dcosQ (1 -x b +Zb-£>) 2 +T} 2 £, 2 

2 



where 



r-NW ,2 e A \V tb \ 2 m t g \Vtb\ 2 m t G 2 F m A w 

C " AW,dl 128(2*)We w ~ NclVudl 4(2^)4 ' (6) 

with 

5 =4, 1= Dt, Z6= 4 (7) 

mf m\v mf 

A convenient way to compute the spin analysing power K p defined in Eq. (jl]) is to evaluate 
the expectation value of cost}: 

(cosr>) = (p-P) = ^JdTcos®= (8) 
We want to compute K p for the following choices of p: 



(i) 


P 


= Pi> 


(ii) 


P 


= Pb, 


(iii) 


P 


= p M 


(iv) 


P 


= h 


(v) 


P 


= T. 



(9) 

In (iv), kj is the direction of the light (non-b-quark) jet with the smallest energy. In the 
leading order calculation of K p , one can simply identify jets with the partons, and thus k, 
denotes the direction of the up-type quark if E u < Ej, and the direction of the down-type 
quark otherwise. Finally, in (v), T denotes the thrust axis [JT7Q. We define the orientation 
of the thrust axis such that T • p# is positive. In leading order, the oriented thrust is given by 
a sign (a • pb), where a denotes the direction of the parton with the largest 3-momentum. 

In Table 1 we list our results. As input we use m t = 175 GeV, nib = 5 GeV, mw = 80.41 
GeV, and Tw = 2.06 GeV. All other constants cancel in the computation of K p . We also 
give numbers for the limiting cases m^ = and Tw — > 0. The latter corresponds to the 
narrow width approximation for the W-boson, i.e. the replacement 

1 71 -h{k 2 -m 2 w ), (10) 



(k 2 -ml/Y + m^T 2 , m w T w 

where k 2 is the squared momentum of the W. One sees that the results are essentially 
insensitive to the bottom quark mass, while keeping the W width changes, e.g., K„ by 
more than 2% as compared to the narrow width approximation. One further comment 
concerning the dependence of K p on the W boson width: In the narrow width approx- 
imation, Tw drops out completely. (This is also true for the QCD-corrected results for 
Kp.) Keeping the Breit-Wigner form for the W propagator, the dependence of K p on the 
W width is extremely small. For example, a 10% change in T\y changes the top quark 
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Table 1 : Born results for spin analysing power of d, b, u , least energetic light jet and 
thrust axis. 





nib = o, rv — > o 


nib — 0, Tw kept 


nib — 5 GeV, TV — ► 


nib = 5 GeV, Tw kept 




1 


1 


1 


1 


J 


-0.40622 


-0.40867 


-0.40553 


-0.40800 


1 


-0.31817 


-0.31091 


-0.31964 


-0.31236 




0.50774 


0.51088 


0.50708 


0.51021 


-0.31712 


-0.31782 


-0.31597 


-0.31671 



hadronic decay rate T by about 10%, but the spin analysing powers are affected only at 
the permill level. Therefore we simply use, both at leading and next-to-leading order, a 
fixed value TV = 2.06 GeV. 

As already mentioned in the introduction, the case of semileptonic decays t — > l + \ib 
follows from the hadronic decay by the identifications l + <-> d , V/ <-> u, and by leaving 
out the factor N c \V ud \ 2 in Eq. @. 



3 QCD corrections 

The computation of the QCD corrections to the decay t(\) — > bdu is a generalization of 
the corresponding computation for f (f ) — > bl + V[ [JT3J3 : For the virtual amplitude, one has 
to add the Wdu gluonic vertex correction (box diagrams do not contribute). The emission 
of real gluons from u and J does not interfere with the gluon emission from t and b and 
are added incoherently. 

We work in d = 4 — 2e space-time dimensions to regularize both soft/collinear and ultra- 
violet singularities. We simplify the virtual amplitude using an anticommuting 75, thus 
respecting the chiral Ward identities. The only divergent part of the amplitude is propor- 
tional to the Born amplitude. The square of the Born amplitude does not depend on d 
if one keeps the W -boson polarization in 4 dimensions. Therefore, we can evaluate all 
necessary traces in 4 dimensions. We also choose to keep the phase space measure dRj, 
in 4 dimensions. Our result for the virtual corrections reads: 

jr virtual _ c a,Q /4V \ £ 1 

dxjdxbd%d cos 6 (1 -x b + Zb-t,) 2 + t\ 2 Z, 2 4n \ mj J T(l-e) 

x [fi{x d -,x b ,Zb) (1 + |P|cos0) +f 2 (xj,x h ,z h )\P\ sin cos x] ,(11) 
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with 



fl(xj,x b ,Zb) 



-2 1 + - ln(co)+ln / (a))+41n(co)ln 



+ 4Li /»( 1 +»-^ 



+ 
+ 



4 



1 +00-X^CO 



S + ln( Zb )+g Wud }x d -(l-x d -z b 



-4Li 2 



1 + co — xt 

1 +CO-XfeCO 



ln(co) 
2z i ln(co 



-Hzb) 



[(l-xj-z b ) 2 +z b (l-z b )] 



■ (2z b -x b + 2xj - xjx b ) + 0(e) 



x b {\ - co) 
1 +co — x^co 



(12) 



f2(x d ,X b ,Z b ) 



^sinQ db x d (l -x d -zb) 



where 



+ 2 Zb 



2(1 -x b +z b ) 
l ~ x d'-Zb 



l-^z b /x 2 



x b 



ln(co) 



ln(co) 



P 

+ 0(e), 



■ln(z&) 



(13) 



(14) 



co 



1-P 
1 + P' 



(15) 



/ 2n e 

8 U 2 



(16) 



and B c j b is the angle between the d and the & in the top quark rest frame. One can easily 
recover the case of semileptonic decays from the above result by leaving out a factor 
Nc\V u d\ 2 and the additional incoherent contribution dl6[ ) from the QCD-corrected Wud 
vertex. Thus we can perform an analytic comparison to the results given in Eqs. (2.7), 



(2.8) of ref. [15]. We find complete agreement by using the well-known correspondence 



1 / 4%fj z 



1117 



1 



r(i 



m I * 

mf 



(17) 



between dimensional regularization in d = 4 — 2e dimensions at a scale /j and regulariza- 
tion by a small gluon mass X used in ref. [[15]]. 

In addition to the virtual corrections we have to include the real corrections which are 
given by the process with an additional gluon: 



t(Pt) -» b(p b ) +u(p u ) + d(pj) + g{p g 



(18) 
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This contribution can also be split into two separate parts, which do not interfere with 
each other: the case where the gluon is emitted from a heavy quark (t or b), and the 
contribution where the gluon is emitted from a light quark. The calculation of the cor- 
responding amplitudes is straightforward. As always in the case of a jet-calculation we 
have to address the question how to cancel the infrared and collinear singularities. For 
the contribution where the gluon is emitted from the secondary fermion line we use the 
dipole formalism [|T5p. Dropping overall factors the subtraction term for this process is 
given by 

with M {pbiPuiPd) being the matrix element at leading order and p, V are defined in 
Eq. (5.3) and Eq. (5.7) of ref. ATS]. The subtraction term integrated over the singular 



region is also given in ref. QT60 and reads 



2% V k2 J r(l-e 



^ + 1 + 10 - n 2 + 0(e) | \M(p b ,p u , Pd )\ 2 . 



(20) 



Comparing with Eq. ([R]) one obtains immediately the cancellation of the soft and collinear 
singularities. For the case where the gluon is emitted from the heavy quark line only a soft 
singularity is present. The collinear singularity is regulated by the finite quark masses. To 
extract the soft singularity we slice the phase space as follows: 

1 = &(2(p r pg) -x min m b m t ) + &(x mia m b m t -2(p r p g )). (21) 

This splits the phase space into a 'resolved' and an 'unresolved' region. The contribution 
from the resolved region is obtained from a numerical integration in 4 dimensions. In the 
unresolved region one can use the soft factorization to approximate the matrix element 
and integrate out the soft gluon. The result is given by: 



1 



dRf\M (p bl p Ul pj,p g )\ 2 ®(x min m b m t -2(pf p g )) 



2m t 

a s C F ( 4k/j 2 x 1 



--ln(zz>) 



2+gln((D) 



2% W^mjJ r(l-£)2m f 
-^(-2p + 2Li 2 (l-«)+ln(«) + lln(w^ (22) 

In order to be consistent with our definition of the phase space measure dR^ for the virtual 
corrections, we have to evaluate the above formula for dR^ =4 . Comparing with Eq. dl2| ) it 
is straightforward to see that the infrared singularities cancel. The numerical implemen- 
tation of the subtraction term given in Eq. ( |T^ ) and the slicing given in Eq. ( pit ) does not 
impose any problem. One should keep in mind that the method presented above allows 
the calculation of arbitrary infrared save observables. 
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4 Numerical results 



While for the three parton final state one can simply identify the partons with the jets, in 
the case of real gluon emission a definition of jets in terms of partons is needed which 
fulfills the condition of infrared-safeness. There are many variations of such definitions, 
and in general the NLO result for K p will depend on the chosen definition. In order to 
study this dependence, we choose two different jet clustering schemes: The E-algorithm 
and the Durham-algorithm. The first step is to compute for all pairs (i, j) of the momenta 
of the final state partons the jet measure yij, which reads 

= E-algorithm, (23) 

mf 

2min{E?,£:?} 

yij = = — — (1 — cos 0(7 ) Durham-algorithm, (24) 

mf 

where 9,y is the angle between parton i and parton j in the top quark rest frame. The 
second step is to recombine the two partons with the smallest yij into a pseudoparticle 
with momentum = pi + pj. Then the direction of the b-jet and the light quark jet with 
the smallest energy can be readily obtained from p^ and the remaining two momenta. For 
bscg events also the direction of the charm jet can be defined if the charm is tagged. The 
construction of the "s-jet" from Meg events is not straightforward: Only events where 
both a &-jet and a c-jet are tagged can be used (in particular, the (rare) events where 
b and c are clustered into a single jet have to be discarded), and for those events the 
remaining third jet is defined to be the "s-jet". This definition includes cases where the s 
is recombined with the b- or c-quark and the remaining third jet consists of a hard gluon 
rather than an s-quark. 

Note that no jet resolution parameter enters in the above definitions; this is not necessary, 
since the leading order process is free from soft and collinear singularities. 

In an experiment which produces top quarks that decay into jets, the first step in the 
analysis is to identify the signal by using a jet finding algorithm and by applying a number 
of cuts. For example, at the Tevatron a cone jet algorithm is used to classify the events 
according to the number of jets. To be more specific, consider top quark pairs where the 
t decays semileptonically and the t decays into jets. Then the event contains at least two 
£-jets and two light jets. For those events with 4 or more jets originating from the t, our 
algorithm should be used in addition to the production-specific jet algorithm, thus leaving 
only events with exactly 3 jets from the t decay. In principle the value for K p can depend 
on the details of the "pre-clustering". This should be studied in Monte-Carlo simulations. 

We can also compute K p for bare b,d (s) and u (c) quarks, since the directions of a 
quark and a quark plus a collinear or soft gluon are identical and thus the condition of 
infrared/collinear safeness is fulfilled. This serves as a benchmark for the realistic case of 
K p for jets. Note however that recombination is needed to define the direction of the least 
energetic light jet. In the case of the thrust axis the concept of a jet is not needed. 
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We write our results to order in the following form 
pO K o _|_ a pi K i pi 

K P = ^ + (XB J r i " =K °P + Wl fO ^ " ^ + ^ = K ' [1 + 6 P CD] + ' (25) 

where denotes the Born result. Table 2 gives our results for K p and §p CD , where we 
use the expanded form of Eq. (|25|). The strong coupling constant is set to a s (m t ) = 0. 108. 



Table 2: QCD-corrected results for spin analysing powers. 





partons 


jets, E-alg. 


jets, D-alg. 


K J 




0.9664(7) 


0.9379(8) 


0.9327(8) 




%] 


-3.36 ±0.07 


-6.21 ±0.08 


-6.73 ±0.08 


K h 




-0.3925(6) 


-0.3907(6) 


-0.3910(6) 


8f° [ 


%] 


-3.80±0.15 


-4.24±0.15 


-4.18±0.15 






-0.3167(6) 


-0.3032(6) 


-0.3054(6) 




%] 


+1.39±0.19 


-2.93±0.19 


-2.22±0.19 


K ; 






0.4736(7) 


0.4734(7) 


8f " [%] 




-7.18±0.13 


-7.21 ±0.13 


Kr 
hf D [ 


%] 


-0.3083(6) 
-2.65 ±0.19 







5 Discussion and conclusions 



Our results listed in Table 2 show that the top-spin analysing powers of the final states in 
non-leptonic top quark decays receive QCD corrections in the range ±1.4% to —7.2%. 
This has to be contrasted with the spin analysing power of the charged lepton in decays 
f(t) -> bl + \f. the QCD corrected result (for m b = 0) reads [|T5p K/ = 1 -0.015a x , i.e. 
the correction is at the permill level. The QCD corrections to the spin analysing power 
of a given quark are much larger due to hard gluon emission from that quark. The spin 
analysing power of jets is smaller than that of the corresponding bare quarks. This effect 
is largest for the "5-jet". We find only a small (at most 0.7%) dependence of the results on 
the jet algorithm. In practice the most important spin analysers are, as far as non-leptonic 
top decays are concerned, the ^-quark jet and the least energetic light (non-b-quark) jet. 
The QCD corrected results are Kb ~ —0.39 and Kj w 0.47. For the £>-jet the difference 
between the parton level result and the jet result is small. The oriented thrust axis, for 
which Kt ~ —0.31, may also serve as a good spin analyser, since it is easily measurable. 
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In summary, we have computed the QCD corrections to the top-spin analysing power of 
jets and the thrust axis in non-leptonic polarized top quark decays. Our results can be 
used in conjunction with the known NLO QCD results for the production of polarized top 
quarks both at lepton and hadron colliders. 
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